Electron localization and entanglement in a two-electron quantum dot 
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Calculations for two electrons in an elliptic quantum dot, using symmetry breaking at the unre- 
stricted Hartree-Fock level and subsequent restoration of the broken parity via projection techniques, 
show that the electrons can localize and form a molecular dimer, described by a Heitler-London- 
type wave function. The calculated singlet-triplet splitting (J) as a function of the magnetic field 
(B) agrees with cotunneling measurements. Knowledge of the dot shape and of J(B) allows de- 
termination of the degree of entanglement in the ground state of the dot, which is of interest for 
the implementation of quantum logic gates. The theoretical value agrees with the experimental 
estimates. 
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Electron localization leading to formation of 
molecular-like structures [so-called Wigner molecules 
(WMs)] within a single circular two-dimensional (2D) 
quantum dot (QD) at zero magnetic field (B) has been 
theoretically predicted to occur Q, Q, as the strength 
of the interaction relative to the zero-point energy 
increases. Of particular interest is a two-electron (2e) 
WM, in light of the proposal [3j for the implementation 
of a solid-state quantum logic gate that employs two 
coupled one-electron QDs (double dot). 

In currently fabricated circular QDs the strength of 
the effective Coulomb repulsion is significantly reduced 
(compared to the value used in Refs. Q, , appropriate 
for bulk GaAs) due to screening from the electrostatic 
ates and the influence of the finite height of the dot 
>]. However, changing the shape of the QD (from a 
circular shape to an elliptical one, and ultimately to a 
quasi-linear one) will reduce the zero-point energy, thus 
enhancing the relative importance of the Coulomb repul- 
sion; this assists in bringing the QD into the strongly cor- 
related regime resulting in electron localization. Here we 
show that this theoretical prediction jfj has indeed been 
recently observed experimentally for two electrons in an 
elliptic lateral QD |7J. We present microscopic calcula- 
tions for two electrons in an elliptic QD specified by the 
parameters of the experimental device [7j . These calcula- 
tions show formation of an electron molecular dimer and 
yield good agreement with the measured J{B) curve (the 
singlet-triplet splitting) when the value of the Coulomb 
repulsion is weakened (by 40%). 

Of special interest for quantum computing is the degree 
of entanglement exhibited by the two-electron molecule 
in its singlet state sB. A measure of entanglement, in- 
troduced in Ref. |8|, is known as the concurrence. To 
date, however, applications H, @ of this measure have 
been restricted to the specific singlet state associated 
with the bonding and antibonding orbitals of a 2e dou- 
ble dot with weak interdot-tunneling coupling. We show 
that our wave- function-based method enables calculation 
of the concurrence of the 2e singlet state in more general 



cases, and in particular for the single elliptic QD of Ref. 
0- This is based on our finding that the concurrence of 
the singlet state is directly related to the degree of parity 
breaking present in the unrestricted Hartree-Fock (UHF) 
orbitals (see below) resulting from our calculations. 

We show that knowledge of the dot shape and the J{B) 
curve allows theoretical determination of the degree of 
entanglement. This supports the experimental assertion 
7] that cotunneling spectroscopy can probe properties of 
the electronic wave function of the QD, and not merely 
its low-energy spectrum. 

The hamiltonian for two 2D interacting electrons is 



H = H(r 1 ) + H(r 2 )+e 2 /( K r 12 ), 



(1) 



where the last term is the Coulomb repulsion, n is the 
dielectric constant, and ri2 = |ri— r2 1 . H(r) is the single- 
particle hamiltonian for an electron in an external per- 
pendicular magnetic field B and an appropriate confine- 
ment potential. For an elliptic QD, the single-particle 
hamiltonian is written as 



H(r) 
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where T = (p - eA/c) 2 /2m* , with A = 0.5(-By, Bx, 0) 
being the vector potential in the symmetric gauge, m* 
is the effective mass and p is the linear momentum of 
the electron. The last term in Eq. J2J is the Zeeman 
interaction with g* being the effective g factor, fis the 
Bohr magneton, and s the spin of an individual electron. 

Our method for solving the two-body problem defined 
by the hamiltoninian consists of two steps. In the first 
step, we solve selfconsistently the ensuing unrestricted 
Hartree-Fock (UHF) equations allowing for breaking of 
the total spin [and parity along the long axis (x-axis) 
of the elliptic QD]. For the S z — solution, this step 
produces two single-electron orbitals ul,r{t) that are lo- 
calized left (L) and right (R) of the center of the QD. At 
this step, the many-body wave function is a single Slater 
determinant ^uhf(1 t,2 j) = \ul(1 T) u -r( 2 I)) made 
out of the two occupied UHF spin-orbitals ml(1 t) = 
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FIG. 1: The singlet-triplet splitting J = E s - E 1 as a func- 
tion of the magnetic field for an elliptic QD with hu x = 1.2 
meV and Tiuj y = 3.3 meV (these values correspond to the 
device of Ref. 0). Lower (solid) curve: The GHL (broken- 
symmetry UHF + restoration of symmetries) result. Upper 
(dashed) curve: The restricted HF (RHF, no symmetry break- 
ing) result. The experimental measurements are denoted by 
the open squares 0. The material parameters used are: 
m*(GaAs)= 0.067m e , re = 22.0, and g* = (see text). The 
calculated values for the concurrence of the singlet GHL state 
(C s , see text) at B = 0, 1.3 T, and 2.5 T are also displayed. In 
the singlet RHF case (upper curve), the concurrence is iden- 
tically zero for all values of B. Note that our sign convention 
for J is opposite to that in Ref. 0. 



UL{vi)a(\) and u R (2 |) = u R (r 2 )P{2), where a((3) de- 
notes the up (down) [j (J,)] spin. This UHF determinant 
is an eigenfunction of the projection S z of the total spin 
S = Si + s 2 , but not of S 2 (or the parity space-reflection 
operator). 

In the second step, we restore the broken parity and 
total-spin symmetries by applying to the UHF determi- 
nant the projection operator [6| P s,t = l^fvji 2 , where the 
operator w\ 2 interchanges the spins of the two electrons; 
the minus sign corresponds to the singlet. The final re- 
sult is a generalized Heitler-London (GHL) two-electron 
wave function \I>Q HL (ri , r 2 ) for the ground-state singlet 
(index s) and first-excited triplet (index t), which uses 
the UHF localized orbitals, 

*GHL( r i: r 2) K (u L (r 1 )u R (r 2 )±u L (r 2 )u R (r 1 ))x s ' t , (3) 

where x s,t = (ct(l)/3(2) =p a(2)/3(l)) is the spin function 
for the 2e singlet and triplet states. 

The use of optimized UHF orbitals in the GHL is suit- 
able for treating single elongated QDs. The GHL is 
equally applicable to double QDs with arbitrary interdot- 
tunncling coupling [6j. In contrast, the Heitler-London 
(HL) treatment [lOj (known also as Valence bond), where 
non-optimized "atomic" orbitals of two isolated QDs are 
used, is appropriate only for the case of a double dot with 
small interdot-tunneling coupling 0, 0] . 

Both the GHL singlet, ^qhl, and the GHL triplet, 
^ghlj cannot be reduced to a single Slater determinant. 



They are always the sum of two Slater determinants, i.e., 
*ghl « M 1 T)«fl(2 I)) T Ml i)u R {2 t)>, andtlms 
they represent states with intrinsic entanglement |3J, li| , 
a property that underlies the operation of the quantum 
logic gate. 

The 2D UHF equations that we are using are described 
in detail in Ref. [1(c)]. For the formulas which enable cal- 
culation of the energies -Eqjj L for the special case of the 
two-electron wave function see Refs. 0,0]. The key 
point is that we exploit an additional variational freedom 
by allowing for different orbitals [here ut,(r) and u R (r)] 
for the two spin directions (up and down). Under ap- 
propriate conditions, the UHF equations do indeed have 
solutions associated with two spatially separated orbitals 
(symmetry breaking is present). We note, however, that 
for sufficiently weak Coulomb repulsion (very large di- 
electric constant k, see below), the two orbitals Ul(t) 
and u R (r) do collapse onto the same single orbital u(r), 
and there is no symmetry breaking. 

Another pertinent point here is that the orbitals 
Ul,r(y) are expanded in a real Cartesian harmonic- 
oscillator basis, namely, 



u l,r( t ) 



k 



L.R 



(4) 



where the index j = (m,n) and = X m (x)Y n (y), 

with X m (Y n ) being the eigenfunctions of the one- 
dimensional oscillator in the x(y) direction with fre- 
quency u! x (ujy). The parity operator V yields VX m {x) — 
(—l) m X m (x), and similarly for Y n (y). The expansion co- 
efficients Cj ' are real for B — and complex for finite 
B. 

We turn now to the interpretation of the measurements 
of Ref. 0- To model the experimental elliptic QD device, 
we take, following Ref. 0, hui x — 1.2 meV and hoJ y = 
3.3 meV. The effective mass of the electron is taken as 
to* = 0.067TO e (GaAs). Since the experiment did not 
resolve the lifting of the triplet degeneracy caused by the 
Zeeman term, we take g* — 0. Using our two step method 
described above, we calculate the singlet-triplet splitting 
Jghl(-B) = Eq HL (B) - Eq HL (B) as a function of the 
magnetic field in the range < B < 2.5 T. Similar to an 
earlier study Q , we find that a weakening of the Coulomb 
repulsion from its value in GaAs (re = 12.9) is required 
in order to reproduce the experimental J(B) curve. The 
effect of screening by the gates can be modeled, to first 
approximation, by increasing re pj.ll3j|. Indeed, with re = 
22.0, very good agreement with the experimental data 
(see Fig. 1) is obtained. In particular, we note the singlet- 
triplet (ST) crossing about 1.3 T, and the backbending 
of the J{B) curve beyond this crossing. 

In Fig. 1, we also plot for re = 22.0 the Jrhf(-S) curve 
[upper (dashed) line] obtained from the restricted HF 
(RHF) calculation, namely for a self-consistent Hartree- 
Fock variation with the restriction that the parity be 
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FIG. 2: Spin-up and spin-down orbitals (modulus square) 
and total electron denities (EDs) associated with the singlet 
state of two electrons in the elliptic QD. The arrows indicate a 
spin-up or a spin-down electron, (a) Doubly occupied orbitals 
of the parity-preserving RHF at B = (left) and B = 2.5 T 
(right), (b) Orbitals of the broken-symmetry UHF at B — 0. 
(c-d) EDs at B = for the UHF (c) and the GHL (d) wave 
functions. The GHL ED is larger in the region between the 
two humps compared to the UHF one; this is indicative of 
covalent bonding. The rest of the parameters are as in Fig. 
1. Lengths in nm and orbital densities in 10 -4 nm~ 2 . 



preserved, so that there is no symmetry breaking and 
both the spin-up and spin-down electrons occupy the 
same spatial orbital [i.e., in this approximation ul(t) = 
Ur(t) — u(r)]. The Jrhf(-B) curve fails substantially 
to reproduce the experimental results. The RHF curve 
disagrees with the experimental trends in two important 
ways: (I) Jrhf(0) is larger than zero, which not only con- 
tradicts the experiment, but also a fundamental theoret- 
ical result [1(b), 14] that states that for two electrons at 
zero magnetic field the singlet is always the ground state, 
and (II) The Jrhf(-B) curve diverges as B increases, un- 
like the Jghl(-B) curve (as well as the experimental ob- 
servation) which bends back after the ST crossing and 
approaches asymptotically the J — line. 

The sharp contrast in the behavior of the RHF and 
GHL results further highlights the significance of the lat- 
ter. Indeed the behavior of the RHF solution reflects the 
independent-particle-model nature of the RHF orbitals 
[see Fig. 2(a)], which leads to a preponderance of the 
exchange contribution, and thus to spontaneous (but er- 
roneous) magnetization of the electron system. On the 
other hand, the GHL results in Fig. 1 reflect the fact that 



symmetry breaking and electron localization [see orbitals 
in Fig. 2(b)] reduce the Coulomb repulsion to a degree 
that compensates for the loss of exchange binding (since 
the localization reduces the orbital overlap) . This is por- 
trayed in formation of a Wigner molecule with the afore- 
mentioned agreement with the experimental J(B) curve. 
The total electron density (ED) of the WM resulting from 
breaking of symmetry [UHF, Fig. 2(c)] and after restora- 
tion of parity (and total-spin) symmetry [GHL, Fig. 2(d)] 
illustrates the electron molecular dimer. We note that 
the UHF solutions at B = exhibit breaking of the par- 
ity for values of re as large as re = 40.0, which indicates 
that the value of re = 22.0 in our calculations places the 
elliptic device of Ref. [j| well within the regime of strong 
electron correlations. 

The asymptotic energetic convergence (beyond the ST 
point) of the singlet and triplet states, i.e., [J(B) — > 
as B — > oo] is a reflection of the dissociation of the 2e 
molecule, since the ground-state energy of two fully spa- 
tially separated electrons (zero overlap) does not depend 
on the total spin. Indeed, orbitals and EDs for finite B 
values are similar to those in Figs. 2(b-d), but with en- 
hanced localization reflected in diminished overlaps with 
increasing B. 

Both the GHL singlet and triplet wave functions [Eq. 
P|)] exhibit entanglement by being the sum of two Slater 
determinants 0, However, unlike the triplet GHL 
state which is maximally entangled (see below), the sin- 
glet GHL wave function may exhibit a reduced degree 
of entanglement, as was found from the Hund-Mulliken 
treatment of the double dot |9j- For purposes of quan- 
tum computing, it is necessary to be able to extract from 
this GHL singlet a quantitative measure of the degree of 
entanglement, e.g., the corresponding value of the con- 
currence C The concurrence, however, was intro- 
duced for the Hilbert space associated with the bonding 
and antibonding orbitals of a double dot molecule [3, li| • 
In our case of a single elliptic QD, the concept of a bond- 
ing and antibonding orbital does not straightforwardly 
apply. Rather, we need to utilize the universal symme- 
try properties of the bonding and antibonding orbitals 
with respect to the parity operator, namely, the bonding 
orbital is symmetric (+), and the antibonding orbital is 
antisymmetric (— ) with respect to reflection about the 
origin of the x-axis. We first notice that, due to the sym- 
metry breaking, the UHF orbitals «z,(r) and ur(t) are 
not eigenstates of the parity operator, and we proceed to 
separate the symmetric $ + (r) and antisymmetric $ - (r) 
components in their expansion given by Eq. Q. That is, 
we write u L ' R (r) oc $ + (r) ±£<I>~ (r); note that $+(r) and 
$~(r) are eigenfunctions of the parity operator. Subse- 
quently, with the use of Eq. (J3J), the GHL singlet can be 
rearranged as follows: 

cx |$+(1 t)$ + (2 I)) - T)#-(2 I)), (5) 

where the coefficient in front of the second determinant, 
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TABLE I: Left: Expansion coefficients Cf' R [See Eq. ®] for 
the broken-symmetry UHF orbitals ul(t) and un(r) at B — 
and k = 22.0. Right: Expansion coefficients for the RHF com- 
mon orbital u(r). The real-Cartesian-harmonic-oscillator ba- 
sis functions are given by <fj (r) = X m (x)Y n (y). The running 
index j stands for the pair (m, n), where m and n denote the 
number of nodes along the x and y directions, respectively. 
A total of K — 79 basis states tpj(r) were used. However, 
only coefficients with absolute value larger than 0.01 are listed 
here. For numbers with a ± in the front, the + corresponds 
to the left (L) orbital and the — to the right (R) one. 



UHF RHF 



3 (m,n) 


°.i 


j (m,n) 


Gj 


1 (0,0) 


0.7990 


1 (0,0) 


0.9710 


2 (1,0) 


±0.5600 


3 (2,0) 


0.2310 


3 (2,0) 


0.2150 


7 (4,0) 


0.0172 


5 (3,0) 


±0.0104 


10 (0,2) 


0.0613 


7 (4,0) 


-0.0180 






10 (0,2) 


0.0253 







rj = £ 2 , is the so-called interaction parameter Know- 
ing T] allows a direct evaluation of the concurrence of the 
singlet state, since C s = 2r]/(l + r] 2 ) 0. Note that $+(r) 
and <f> _ (r) are properly normalized and that they are by 
construction orthogonal. 

For the GHL triplet, one obtains an expression inde- 
pendent of the interaction parameter 77, i.e., 

*GHL « T)*"(2 I)) + |$ + (1 i)*~(2 T)>, (6) 

which is a maximally (C* = 1) entangled state. Note 
that underlying the analysis of Ref. 7] is a conjecture 
that wave functions of the form given in Eqs. 10) and @ 
describe the two electrons in the elliptic QD. 

To make things more concrete, we display in TABLE 
I for B = (and k = 22.0) the coefficients C^' R [see 
Eq. 0] that specify the broken-symmetry UHF orbitals 
in the real Cartesian harmonic-oscillator basis. Indeed, 
we find that there are contributions from both sym- 
metric and antisymmetric basis functions X m (x)Y n (y) 
along the x-axis, since both even and odd m indices 
are present. The n indices are all even, since there 
is no symmetry breaking along the y-axis. Naturally, 
all the m indices in the expansion of the RHF com- 
mon u(r) orbital are even, since in this case the par- 
ity is preserved. One can immediately check that £ — 

M \r L \ 2 /v x \n L \ 2 Y 

(m odd) I j I / Z—/j (m even) I j I J 

For the RHF singlet £ = 0, since there no coefficients 
with odd m indices (TABLE I). For the GHL singlet, 
we calculate the interaction parameter rj = £ 2 and the 
concurrence C s for the device of Ref. in the magnetic- 
field range < B < 2.5 T. We find that starting with 
rj = 0.46 (C s = 0.76) at B = 0, the interaction param- 
eter (singlet-state concurrence) increases monotonically 
to 7] = 0.65 (C s = 0.92) at B = 2.5 T. At the interme- 
diate value corresponding to the ST transition (B = 1.3 



T), we find rj = 0.54 (C s = 0.83). Our B = theoretical 
result for rj and C s are in remarkable agreement with the 
experimental estimates of r/ — 0.5 ± 0.1 and C s = 0.8, 
which were based solely on conductance measurements 
below the ST transition (i.e., near B = 0). Our theoret- 
ical values for the singlet-state concurrence are marked 
on the GHL curve in Fig. 1. We also display in this figure 
the value of C s for the RHF single-determinant solution 
(upper curve). The RHF value for this quantity (and 
also for rj) vanishes identically for all B values. Thus the 
shape (and values) of the experimental J(B) curve por- 
tray directly the degree of electron localization and the 
entanglement associated with the singlet state in Eq. 
Note that our calculated results are based on the ener- 
getics of the singlet-triplet splitting alone, and thus they 
provide a reliable and independent alternative method 
for extracting the degree of singlet-state entanglement 
for all values of B. 

In conclusion, we have shown formation of an electron 
molecular dimer in an elliptic QD (Fig. 2) for screened 
intcrclcctron repulsion characterized by a singlet-triplet 
splitting J(B) that agrees with experiment (Fig. 1). Fur- 
thermore we showed that, from a knowledge of the dot 
shape and of J{B), theoretical analysis along the lines 
introduced here allows probing of the correlated ground- 
state wave function and determination of its degree of 
entanglement. Such information is of great value for the 
implementation of solid-state quantum logic gates. It is 
also of interest to quantum information theory in general 
[I5I Il6| . The present theoretical method and analysis 
can be straightforwardly extended to double dots with 
arbitrary interdot-tunneling coupling. 
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